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Abstract: 

We discuss the exact renormalization group or flow equation for the effective action 
and its decomposition into one particle irreducible N point functions. With the 
help of a truncated flow equation for the four point function we study the bound 
state problem for scalar fields. A combination of analytic and numerical methods 
is proposed, which is applied to the Wick-Cutkosky model and a QCD-motivated 
interaction. We present results for the bound state masses and the Bethe-Salpeter 
wave function. 
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1 Introduction 



For many considerations in quantum field theory it is a helpful idea to integrate out 
high frequency or short distance modes, and to study a scale-dependent effective 
theory for the remaining low frequency modes. This concept was first made con- 
crete on the lattice within the framework of block spin transformations |l| . When 
applied to continuum quantum field theory, it can be cast into the form of exact 
renormalization group equations p[-|10[, which describe the cutoff dependence of 
the effective theory in a compact way. Polchinski and others 0-0 employed these 
ideas in order to study the dependence of physical Green functions on the UV cutoff 
with the aim to simplify perturbative proofs of renormalizability. The integration 
of the corresponding flow equations with respect to an infrared cutoff, on the other 
hand, provides a new exact method to compute effective low energy Lagrangians 
T0[] or average actions [|Il] . Within perturbative or 1 /N expansions this scheme has 
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been applied successfully to the computation of effective potentials and "running 
wave function normalizations and allows, e.g., to investigate cricital phenomena in 



two and three dimensions |TT[|, [T^, high temperature phase transitions |T3[ and 
universality within the Higgs top system [jlO| . 

These results employed, to a large extent, an expansion of the effective action in 
powers of momenta (as it is appropriate for investigations of the effective potential). 
The first purpose of the present paper is the derivation of an expansion of the exact 
flow equations for the effective action in powers of fields, keeping the momentum 
dependence exact. This way one finds an infinite set of coupled differential equation 
for one-particle irreducible N point functions. 

Then we concentrate on the four point function, whose singularities as a function 
of the cm. energy squared s encode informations on possible bound states of theory. 
We study a truncation of the infinite set of flow equations, which corresponds to 
the ladder approximation. It leads to a simple flow equation including just the four 
point function itself. Already in this approximation we see, that the flow equation 
allows to find both the masses and the wave functions of bound states in a theory, 
and the equivalence to the Bethe-Salpeter equation can be shown. 

Next we discuss, as an essential step towards numerical integrations of the flow 
equations, the Laplace transform of the four point function with respect to the 
Lorentz invariant products of momentum variables. Numerical methods can then 
be introduced after discretization in this "Laplace space". Integration of the flow 
equation becomes a simple algorithm to update the corresponding "Laplace lattice" 
consecutively. The informations on bound states can be obtained from the numerical 
results in combination with analytic methods. 

As a first application we study the Wick-Cutkosky model, which contains two 
massive complex scalars interacting through the exchange of a massless real scalar 
and ressembles, in the nonrelativistic limit, to the Positronium problem. We find 
good agreement between our results, both on the couphng constant dependence of 
the bound state mass and on the bound state wave function, and known formu- 
las. Our method does not make use, however, of the 0(4) symmetry of the Wick- 
Cutkosky model and can straightforwardly be applied to the bound state problem in 
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the case of arbitrary interactions. We present formulas, which allow to study inter- 
actions which correspond to potentials, in the non-relativistic limit, of the form 
with a arbitrary. We end up with results for a QCD-motivated potential including 
a linearly rising confining part. 



2 Flow Equations for N Point Functions 

We will present the general features of the flow equations in the context of a single 
scalar field for simplicity, and we work in Euclidean space. The aim of the flow 
equations is the computation of Green functions within a theory, which is regularized 
in the UV in terms of a cutoff A in the propagator. In addition an infrared cutoff k is 
introduced below, and an important role is played by the corresponding propagator 

pt{^') - {R^M'))-' = (2.1) 

_|_ ^/ 

with 

hk{q^) 1 for < ^ 

hk{q^) for q^ > A;^ (2.2) 

The starting point is the UV-regularized generating functional of connected 
Green functions G^{J), which can be represented as 

g-G'^(J) = X [ V(l)e-^^'^'^^'t>)-stM)+iJ,'t>) _ (2.3) 



Here (J, 0) etc. is a short-hand notation for 

and we have represented the kinetic term in terms of the inverse UV-regularized 
propagator i?g (g^). An alternative and useful representation of G^{J) is given by 



with 



[2.5) 



Next we introduce an infrared cutoff k and define an effective interaction Sint{4>, k) 

by 

where is given as in (|2.6|) with Pq replace by P^. Clearly Sint{4>, k) is equal to 
the bare interaction S^^^^cp) at k = A, and from ( p.5|) we have 
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The flow equation for Sint{4>,k) follows easily after differentiation of the defining 
equation ( |2.7| ) with respect to k (subsequently we us dk = d/dk"^) 0, [|T0| : 

dkSU<P,k) --J j^APk iq ) • - 50F^/- ^^-^^ 

The generating functional of connected Green functions including the infrared 
cutoff k, G^{J), is related to Sint{(p, k) through 

Sint{m\^^p.,j = Gt{J) + \{J,PtJ)- (2.10) 

Inserting (p.lO|) into (|2.9|) and taking into account that J is implicitly fc-dependent 
because of the = ■ J prescription in (p.lCI| ), one finds that G^{J) satisfies a 
similar flow equation: 

d.G, (J) = -- J j^ARk ■ |^7(^)^7(Z^ - 7j(^5jF^/- ^^-^^^ 

In ( p.ll|) we have neglected a J-independent term in G^{J). For many considerations 
in quantum field theory, however, it is more convenient to deal with the generating 
functional of one-particle irreducible Green functions T (</?). Its constant part gives 
the effective potential, and also all other informations on a theory are encoded in 
r(<^) in a more compact way. In the presence of a UV cutoff A and an infrared 
cutoff k the effective action T^{ip) is given by the Legendre transform of G^{J), 

Tt{^) = GUJ) + {J,^). (2.12) 

Inserting the Legendre transformation into ( p.ll| ) it is possible to derive the flow 
equation for T^{ip). Thereby the implicit /^-dependences of J and ip have to be taken 
into account, but at the end they cancel and one is left with the simple result 

Such an equation has also been found be Wetterich and, in a different context 
(as a differential equation with respect to the UV cutoff A instead of k) in 0. If 
one splits off a bare kinetic part of T^{ip), 

n{v) = l{v,Rt^) + Tt{^), (2.14) 



one obtains a flow equation for f ^ (v?) of the form 

*f '(^' - \ I ('^) (''^ ('^) + " • P-^^' 

In general, however, will still contain terms quadratic in ip. The boundary 
condition for /c — > A is most easily given for Tk{'p>)'- From the relation between 
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T^{(p) and Sint{4>,k) via eqs. ( ^.121 ) and ( |2.1(J| ), and a careful consideration of the 



limit for k ^ A, one finds 



f i(^) = StM ^ (2.16) 



where Si^^{(j)) is the bare action of eqs. ( |2.3D or 

Let us now consider an expansion of in powers of with one particle irre- 
ducible point functions r^r as coefficients. 

For simplicity we will assume a discrete symmetry —ip such that only even 
powers of ip appear. Instead of writing the convolutions over momenta, we switch 
to an tensor notation (as appropriate in the case of discretized momenta as, e.g., on 
a torus) according to 

^(g) ^ ^„ r^(gi...g^) ^ F- '^-,/ ^^T. ,etc. (2.17) 

Thus we have 

1 • • 1 

+ —^l^iVjVkVlVmVn^i + ••• (2.18) 

Inserting this expansion into the flow equation ( |2.13| ) and ordering the result ac- 
cording to powers of y?, one obtains the following inflnite system of equations: 



2 ~ ztm~ 'i ~ 2nh\ 



1. 



op— 1 -popij-p — \ -pqrhl-p — l P^l -popijhl-p~\_ 

2no'' 4 2pq'- 4 2rm ~ 2 "^"-o ^ 2pm 



(2.19) 

The r.h.s. of the flow equations for the individual N point functions Fat have the 
following properties: In contrast to standard (3 functions each of them is exact, there 
are no higher order terms beyond the ones shown explicitly. They have a simple 
diagrammatic interpretation shown in flg. 1. The flow equation for F^v includes 
always a term involving F7V+2 on its r.h.s, thus a truncation of the series leaves no 
exact result. On the other hand, the iterative solution of these differential equations 
generates the loop expansion of standard perturbation theory, because the r.h.s. of 
each equation involves effectively an overall factor of h. 



3 Flow Equation for the Four Point Function 

In this section we will consider a truncation of the flow equation of the four point 
function, which corresponds to the ladder approximation in the Bethe-Salpeter 
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framework. First, we will neglect the contributions to the flow of the full inverse 
propagator r2 induced by the interactions and identify r2 with the bare (regularized) 
inverse propag ator, = (or T^^ = P^). Second, we neglect the contribution 
of Tq to the flow of which amount to taking only bubble-type diagrams into 
account. 

We will slightly change the field content, namely include two distinct complex 
scalar fields ^Jq, ipi, with identical mass m for simplicity. The expansion of the effec- 
tive action reads 

r^(^„ = ^iRf ^cpi + ^IX' M + + - (3.1) 

and the flow equation for r4 becomes in the present approximation 

o -pih _ r,a pAm -pnh pAo -pip ('i 0\ 

(Jki- 4jl — ^(-'fc-' fc n'- Aol^k p'- ijm- 

With momenta written explicitly this equation reads 

<9fcr4(pi,P2,P3,P4) = -2 j -^;^T4{puP2,q,-q-pi-p2) 

■r,{-q, q + pi+ P2, P3, P4) ■ dkPtiq') ■ Pkiiq + Pi + P2?). (3.3) 

The equivalence to the Bethe-Salpeter equation can be established as follows: 
Let us introduce a two-particle propagator P^^) through 

p(2)ij _ pAi pAj /o A\ 

mn — fc m-* fc rf l-J'^J 

Then eq. ( p.2|) can formally be written as 

dkV^ = -14 ® 9fcP(') ® 14 (3.5) 

with the formal solution 

pA 

= 1 + Pii ® n ^^'^^ 

where r4 is the boundary value or bare four point function. In our convention, the 
non-amputated four point function r4, which also includes a disconnected piece, is 
related to r4 through 

f 4 = -p(2) ® 14 ® P(') + P(') (3.7) 
and satisfies the Bethe-Salpeter equation [|1^ 

f 4 = P(2) + p(2) ® K ® f 4 (3.8) 
where K denotes the interaction kernel. The formal solution of (B.SI) reads 



P<^. ^ ,3.9) 



i-/s:®p(2) i-/s:®p{2) ■ 

Using eq. (|3.7|) , one finds that this solution of the B.-S. equation coincides with the 
solution ( p. 61) of the flow equation after identifying K with — r4 . Thus, although 
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these equations are quite different, they have the same physical content. We expect 



that a general equivalence between the full set of flow equations ( p.l9|) , beyond 
the ladder approximation, and the full set of Schwinger- Dyson equations can be 
established. 

This discussion showed already an equivalence between the boundary value of 
the running four point function r4 at A; = A, denoted by r4 in eq. (^.61), and a 
general interaction kernel K. In fact, r4 does not necessarily have to be a bare 
coupling of a theory with only ipa, V^fe fields, but could be the result of interactions 
with different fields which have already been completely integrated out. An example 



is given by the Wick-Cutkosky model 14 1, which starts out with an additional real 



massless scalar field (j) beyond the fields y^^, v?6- The original interaction involves no 
bare four point coupling as in eq. (|3.1| ), but just trilinear interactions of the form 



int 



9<P{fifa + fyb)- (3.10) 



It is straightforward to integrate out the field 0, which gives rise to an effective 
(nonlocal) four point function for the field (pa, fb as in (|3.1|) with 



r4(Pl,P2,P3,P4) = 7 ^ = 7 T^- (3-11) 

The dynamics of the remaining fields ipa-, (fib in the presence of the interaction ( |3.11|) 
can now either be obtained by means of the B.S. equation or, as we propose here, by 
integrating the flow equations for r4 as given by eq. ( p.3|) with (|3.11| ) as boundary 



condition. It is easy to imagine more general boundary conditions obtained, e.g., 
after integrating out photons, gluons, or the hke. 

In order to deal with eq. (|3.3|) we next observe that, since the framework is 
completely Lorentz-covariant, r4 should be a function of Lorentz- invariant products 
of momenta only. Off-shell six independent invariants can be formed, which we 
denote by 

S= {Pl+ P2Y = {P3 + Pif, t= {P2- Pzf = {P2 - Pif 

Vi=pI+ pI, V2=pI+ pI, Wi=pI- pI, w2=pI- pI- (3.12) 

Although r4 will in general be an arbitrary function of those six variables, we will 
first study the case where r4 can be written as a product as follows: 

14 = /(f 1, wi) ■ Dk{s) ■ f{v2, W2). (3.13) 



If one inserts this ansatz into the flow equation (|3.3| ), one finds that it is stable 
and thus corresponds to a "fix point" . Furthermore variables can be separated, and 
one finds that the function / can depend on k only via an overall factor, which can 
be included in Dk{s) by definition. This latter function satisfies a flow equation of 
the form 

dkDki-s) = -2Dl{s) ■ FkisJ) (3.14) 
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with 



(3.15) 

The general solution of the flow equation ( ^.141 ) is easily written down in terms of 



an arbitrary boundary condition -Dfc'(s) at k = k': 

Dkis) = ^4^^ (3-16) 

On the one hand, this shows that any problem where r4 starts out to be of the 
form (|3.13|) is analytically solvable. On the other hand, iff the theory contains a 
bound state in the channel {ipa,(pb) — > (v^a,V^fe), which manifests itself as a pole in 
s of we actually expect that r4 approaches the form ( p.l3| ) for A; — > and s in 
the vincinity of the pole. Then the function / corresponds to the amputated Bethe- 
Salpeter wave function of the bound state, whereas -Dfc(s) for k ^ corresponds to 
the bound state propagator and contains the information on the position of the pole. 
In fact we will see that the factorized form ( p.l3|) of r4 is an infrared-attractive "fix 



point" of the fiow equations, which will be of great help in extracting informations 
on possible bound states analytically. 

In order to proceed towards concrete calculations we have to make a choice on 
the form of the regularization of the propagator or on the function hk{q^) of eqs. 
(p.l|) and (|2.2|). The following choice turns out to be particularly useful: 



or 



hk{q') = e~^^ (3.17) 



y -TMi. y -rill. 

Ptd) = ^-^-^-^ ^ (3.18) 



q2 _|_ ^2 Jl/A2 

This regularized propagator has the following properties: For general nonvanishing 
fc^, there is no pole at q^ = — m^, instead it behaves like l/k"^, at q^ = — m^, for 
k"^ ^ 0. Beyond the "pole", for < — m^, it increases exponentially for k"^ 0. 
The absence of singularities at finite k"^ is also clear from the boundedness of the 
Feynman parameter integral da. 

Similarly we have to make a choice on the parametrization of the dependence 
of r4 on the six Lorentz invariants shown in ( ^.12| ). First we make a technical 
simplification, namely we neglect the dependence of r4 on the variables Wi,W2 of 
( ^.12| ). In the Wick-Cutkosky model r4 does not depend on Wi,W2 neither at its 
starting point at A; — A, given by eq. ( p.ll| ) in the form g'^/t, nor for ^ in the 



weak coupling limit. In this limit the wave function / is known to depend on v only. 

Now it turns out to be very convenient to switch from r4(s, t, vi, V2) to its Laplace 
transform with respect to the three variables t,Vi and V2'. 

POO 

Ti{s,t,vi,V2) = dlodkdl2Ck{s,lo,h,l2)e-'''-'''"-'''" (3.19) 
Jo 
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After inserting this Laplace transformation together with the parametrization 
of the propagator into the flow equation (pl3|), one finds that the d'^q integration can 
easily be performed, and one is left with the following flow equation for Ck'- 

dkCkis,lo,h,h) = — / d3 dnodnidn2dmodmidm2 
8n^k^ Jk^/h? Jo 

B B IB 

■Sih TTT^ "^2) ■ Ck{s,no,ni,n2) ■ Cfc(s, mo, mi, ma) (3.20) 

ZJd 

with 

6 = 1 + P(n2 + mi), h' = 13 + k'^{n2 + m^, B = b + b' + P{no + mo). (3.21) 

Of course the previously discussed simplification in the case of factorization remain 
valid: If r4 factorized as in eq. ( p.l3|) , Ck factorizes as 

Ck{s, lo, h, I2) = 6{lo)f{h)Dk{s)f{l2) (3.22) 

where /(/) is the Laplace transform of the wave function f{v). The solution ( p.l6| ) 
for Dk{s) remains the same with Fk{s, f) now given by 

^'=(^'/) = T^i/,.^/5X (3.23) 

and 6 = 1 + P{h + k), b' = (3 + P{h + h), B = b + b'. The right-hand side of 
eq. ( ^.2(J ) shows already an expected analytic property of the four point function: 
If we continue the variable s into the Minkowskian regime (s < 0), we expect a cut 
at s = — 4m^. Indeed, in the limit k'^ ^ and for /3 — > 1 the exponent on the r.h.s. 



of eq. ( |3.20| ) becomes — (2m^ + s/2)/k'^ and thus explodes for s < —4m 
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In the case of general interactions, i.e. general boundary conditions for r4 or 
Ck at k = A, the flow equations either in the form ( |3.3| ) or in the form ( |3.20| ) 
cannot be solved analytically. Now we propose to solve eq. (p.20|) numerically 
after discretization of the variables lo,h and k- Ck{s,lo,h,h) becomes a function 
Ck{s)iQij^i2, which lives on a three-dimensional lattice, where the integers io,ii,i2 
denote the lattice points. After discretization of the variable k^ eq. ( ^.2(J| ) provides 
an algorithm on how to update the lattice in each step from Ck{s) to Ck-Ak{s). 

It is sensible to start this algorithm with values of k large compared with the 
scale m of the problem; choosing the starting point A to be A^ ~ lOm^ turns out 
to be sufficient for our purposes. At the beginning one has to fix Ca, the Laplace 
transform of r4 , and the variable s. 

In the presence of a bound state with mass M we expect the following behaviour 
of Ck in analogy to the behavior of the regularized propag ator Pfc^ of eq. j ^JEj ): 
For s > — Ck remains finite for k ^ 0, whereas for s ^ — Ck approaches 
a factorized form as in eq. ( |3.22[ ), and the factor Dk{s) diverges for s < — and 



k ^ 0. A divergence of Dk for — thus indicates, that the variable s satisfies 
s < — M^. This behaviour is indeed the result of our numerical investigations. 
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Since the factorized form ( |3.22| ) turns out to be infrared stable, eq. ( p. 161 ) can 



help us to decide whether, at a given value of s, Ck^o diverges or not: If we observe 



that for a certain value of A; = k', Ck' assumes the form (|3.22| ), we can read off the 
function /(/) and the value of D^/^s). Then we can compute the function F^i^s, f) 
of eq. ( |3.23| ), insert it into eq. ( p.l6| ), and we see immediately that -D^'^g diverges if 

2Dkis) r dPF^,{s,f) > 1. (3.24) 
Jo 

There is thus no need to integrate numerically down to A; = 0, but the following 
procedure turns out to be sufficient: We integrate numerically the flow equation for 
Cfc, at a given value of s, down to k'^ = k'"^ ~ m^/10. Then we check whether Ck has 
assumed the factorized form ( p.22| ). If not, this remains so even for smaller values 
of k"^ (as we have tested numerically) and we conclude, that the corresponding value 
of s satisfies s > — M^, where M is the mass of the lowest lying bound state. Thus 
we decrease s and start again, until we get to a value of s where Ck'{s) factorizes. 
Then we read off Dk'{s) and /(/) and check with the help of eq. ( p.24|) whether 
Dk^o{s) diverges, i.e. whether s satisfies s < M^. By repeating this procedure for 
different values of s, we can pin down as accurately as we like, and of course we 
have also obtained the Laplace-transformed wave function f{l). 

The general picture of the use of ffow equations is the following: For k"^ ^ m^, 
the r.h.s. of the flow equations is small except for Green functions which correspond 
to marginal or relevant operators, which are thus the only ones to vary signiflcantly 
(but which are not present in the present approximation). In this regime, at least 
for small couplings, the result of the integration of the flow equaitons could also 
be obtained by standard perturbative methods, eventually improved by the use 
of any mass- independent standard renormalization group. For k"^ ~ m^, however, 
dynamical effects like the formation of bound states take place, which generally 
require the use of numerical methods. For k"^ -C m^, the r.h. sides of the flow 
equations are again small (damped exponentially) except for external momenta, 
where Green functions develop singularities like poles or cuts. This is a feature 



both of the bare propagator (|3.18|) and of the four point function discussed here. 
The development of singularities, however, has simple universal properties like the 
factorized behaviour ( |3.13| ) of the four point function, which allow it to be treated 
analytically again. 

The general procedure of combined numerical and analytic methods for the 
search for bound states as described above will be applied to some concrete models 
in the next section. 



4 Results for Bound States 

In this chapter we will study the bound state problem in the case of two complex 
scalar flelds using the methods derived in section 3. We integrate the flow equations 
for the four point function r4 in an approximation, where the scalar propagator 
remains the free one, and where the contribution of the six point function to the 
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r.h.s. of the flow equation for r4 is neglected. After restricting the momentum 



dependence of r4 to the four invariants s, t, vi and V2 of (|3.12| ) and after the Laplace 
transformation with respect to t,vi and V2, the flow equation has the form (|3.20| ). 

The information on the underlying model, i.e. the bare interaction r4 or the 
kernel of the corresponding Bethe-Salpeter equation, is entirely and only specified 
by the boundary condition at k = A. In the case of the Wick-Cutkosky model, 
which we study first, this boundary condition is given by r4 = g'^/t (see eq. ( |3.11|) 
or, in terms of the Laplace transform Ck of ( |3.19| ), by 



Cj,{s,loJiJ2)=g'S{h)6{k). (4.1) 

As described in the previous chapter, we discretisize the variables Iq, li and I2. We 
used lattice sizes up to 20 in all three directions. We integrate the flow equation 
( p.20| ) numerically from /c^ = A^ = lOm^ down to k"^ = k''^ = m^/10 at fixed cm. 
energy s with —Aw? < s < 0. At /c^ = A;'^ we check, whether Ck' has at least 



approximately assumed the factorized form (|3.22| ). If this is the case, we extract the 



function /(/) and the value of Dk'{s). Then we compute the l.h.s. side of ( ^.24| ) 



still at fixed s, in order to check whether s satisfies s > — or s < — M^, where 
M is the mass of the lowest-lying bound state. 

Indeed this method works very well in the present model. If s is chosen close to 
or below — M^, the numerical deviations of Ck, from the factorized form (|3.22|) , are 
less than one percent. For such values of s, Ck is seen to increase rapidly towards 
small values of k. For larger values of s, on the other hand, Ck is seen to remain 
finite for k ^ 0. If we would take Ck'{s, 0, 0, 0) > 10^ as a criterium for s < — M^, 
this would actually match our criterium for s < — based on eq. (|3.24|) already 
within a few percent of the value of s. 

Within the Wick-Cutkosky model it is convenient to study the dependence of 
the mass M of the lowest lying bound state on a coupling A defined by 

A = g^/lGirm^. (4.2) 

Analytic results are known in the weak coupling limit A — >^ 0, where 

~ (4 - A2)m^ (4.3) 



and in the case of a vanishing bound state mass [14|: 



M = for A = 271. (4.4) 



A result of a numerical solution of the Bethe-Salpeter equation can be found in |15|: 



M = (2 - .082)m for A = 1. (4.5) 

In fig. 2 we plot our results for in units of versus A. We indicate, for A < .3, 
eq. ( |4.3| ) as a short line, and the two results ( [4.4| ) and ( [4.5|) as crosses. We show 
our results as error bars, which were obtained on a 20^ lattice. The error bars are 
due to varying the prescriptions on how to discretisize the variables Iq, h and I2, and 
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due to changing the physical size of the lattice: The ranges of Iq, h and I2 are given 
by < /o < lomax and < li,l2 < ^max, and we varied lomax between 10 and 40, 
and Imax between 6 and 10 (in units of m~^). (Larger values of Imax are irrelevant, 
because Cfc(s, /q, /i, decreases exponentially with li and I2, see below.) 

In the case of large values of the coupling A it should be remembered that we 
neglected the dependence of r4 on the momentum variables Wi , W2 of eq. ( |3.12D . This 



is known to be a good approximation, for r4 with k <ti A, only in the nonrelativistic 
or weak coupling limit. Therefore the deviations of our results from the known 
behaviour ([4.4|) for A — ^> 27r (strong coupling) are understandable. 

In the extreme weak coupling limit A ^ 0, on the other hand, the physical size 
of a bound state, in ordinary space, is known to increase. In terms of our "Laplace 
space", spanned by the variables /j, large distance also correspond to large values of 
li. Since this "Laplace space" necessarily has a finite size, it is again understandable 
that this method cannot describe phenomena accurately, where very large distances 
play an important role, as in the case of very weakly coupled bound states. 

For moderate values of A we find a large region, however, where our results seem 
to be free of those systematic errors and match the known ones suprising well. The 
dependence of our results on the lattice size might be of some interest. Therefore 
we plot, for A = 1, our results for versus the lattice size in fig. 3. (The error 
bars are of the same origin as discussed above.) We see the steady convergence, 
with increasing lattice size, towards the known result, which makes us believe that 
the method has no systematic limitations beyond the ones mentioned before. 

The procedure required already to extract the function /(/) out of Ck', which 
satisfies ( p.22| ) for s near — M^. Up to a normalization this function is the Laplace 
transform of the amputated Bethe-Salpeter wave function. In the weak coupling 
limit of the Wick-Cutkosky model this amputated wave function is known to be 



JXpIpD = 1/{pI+pI + 2^2) = l/{v, + 2m2) (4.6) 

with Vi as in ( p.l2| ). Accordingly /(/) should read 

/(/) =exp(-2m2/). (4.7) 

In fig. 4 we plot our result for log [/(/)], at A = 1, versus / in units of for a 
lattice size of 10. We normalized our result with respect to the known one at / = 1, 
and we see very nicely the agreement with the known exponential decrease, until 
finite lattice size effects start to play a role. 

It should be noted that, in order to obtain all the results of figs. 2 to 4 from the 
numerical integration of the flow equation (|3.20|) , a computer time on a work station 
of only a few minutes was needed. Furthermore no use was made of the 0(4) symme- 



try of the Wick-Cutkosky model |14], which was exploited in the numerical solution 



of the Bethe-Salpeter equation in |T5[. Thus the method can straightforwardly be 
applied to the relativistic bound state problem with arbitrary interactions. 

Of particular interest are interaction kernels, which are motivated by QCD and 
which correspond to a confining potential in the nonrelativistic limit. Let us first 
discuss potentials V{r), which behave like r" as a function of the distance r in 
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ordinary three-dimensional space. The corresponding Lorentz-invariant momentum- 
dependent interaction kernel is actually not unique; it becomes unique, however, if 
one assumes that it depends only on the momentum transfer = t. We adopt the 
notation to our procedure, where the interaction kernel coincides with the boundary 
value of the four point function; then the relation between V{r) and T^{t) reads 

r^(t) = 723+V/V3-"^ffl| for V(r)=7r". (4.8) 

In our formalism we actually need the Laplace transform C\ in analogy to ( |4.1| ); its 
corresponding form is given by Ca(s, Iq, li, I2) = C\{lo)6{li)6{l2) with 

93+0^3/2 

Ca{Io) = l^^^lo' ■ (4.9) 

The possible singularities in a in eqs. ( [4.8| ), ( [4.9| ) originate from analytic continua- 
tion in a, which is required if the functions V{r) or Ca(/o) are not well behaved at 
large distances r — cxd or Zq — ^ cxd. Within our procedure, however, Ca(/o) will only 
be nonvanishing for a finite range of Iq due to the finite size of the lattice. From the 
general relation between V{r) and Ca(/o), 

fM = ii;r§C^A(We- w.th ,„ = ^. (4.10) 



one finds the following: If Ca(/o) vanishes for Iq > a, V{r) decays exponentially for 
S> a. On the other hand, for <^ a, V{r) and C\{Iq) are related as given by eqs. 
([4.8[), ( [4.9| ) up to a constant in V{r) depending on a. Such a constant is actually 
quite welcome; mesonic spectra in QCD, in the nonrelativistic limit, are typically 
derived from potentials of the form 

V{r) = \r- ^aQCDT-' - c (4.11) 

with a constant c of 0(1 GeV) p6[| . 

We applied our method to relativistic interactions, which are motivated by the 
potential ([4.11| ). Of course we are dealing with scalars instead of fermionic quarks, 
but scalars are known to reproduce the spectrum of fermionic bound states in the 
nonrelativistic limit (where spin effects become negligible), provided the Coulomb 
interaction induced by vector boson exchange is multiplied by a factor 4. 

Thus we integrated the flow equations ( ^.21 ) with a boundary condition Ca(/o) of 
the form 

Trr^ 4 
— -/o)(87rA/o + 167r--< 

which reproduces ( [4.11|) for r <^ c/A according to (|4.10| ). For the string tension A 



Ca(/o) = -0{^ - /o)(87rA/o + 167r ■ -aQco), (4.12) 



we choose .25 GeV^, .3 for agcD and 1 GeV for the constant c. In table 1 we give 
our results for the bound state mass M for different values of the scalar mass m in 
the range 1.5 GeV to 5.5 GeV. Again the results agree with general expectations: 
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Small, but increasing binding energy towards larger scalar masses, and vanishing 
binding energy for m ^ IGeV (Our present method does not allow us to find 
bound states above threshold or with M > 2m, since for s < —Am? the behaviour 
of the four point function is dominated by the two-particle cut, which reveals itself 
as a singularity for — 0.) The increasing errors towards larger scalar masses and 
binding energies likely indicate that here the neglect of the dependence of r4 on the 
momentum variables Wi,W2 is not a good approximation. Nevertheless the results 
show certainly the feasibility of the method in the case of general interactions. 



5 Conclusions 

In this paper we presented a derivation of the exact flow equation ( p.l3| ), which 
allows to compute effective low energy actions in terms of arbitrary high energy 
actions. It corresponds to an infinite set of fiow equations for one particle irreducible 
Green functions, where each equation is again exact. Solutions to a finite subset of 
these equations are already nonperturbative in and correspond to summations of 
certain sets of Feynman diagrams. They contain much more information, however, 
than the finite number of running coupling constants in standard renormalization 
theory. 

Whereas the practical use of fiow equations, often in different formulations, has 
been demonstrated before in context of expansions of effective actions in powers of 
momenta |l^-|[T^, we concentrated here on the bound state problem. Our tool was a 
truncated flow equation for the four point function, corresponding to the summation 
of "ladder type" diagrams, with the full dependence on the momenta left intact. 

In principle, the effective low energy action and hence the low energy four point 
function contains all informations available by summing perturbation theory. We 
presented some technical tools, which allow to make practical use of this principle. 
These tools involved the parametrization of the momentum dependence of Green 
function in terms of its Laplace transform with respect to the Lorentz-invariant 
products of momenta (except for the cm. energy s). This procedure is obviously 
manifest Lorentz-invariant, and simplifles the r.h.s. of the flow equations such that 
they become accessible to a numerical treatment. 

One has not to rely completely on numerical methods, however: In the presence 
of nontrivial analytic structures, such as bound state poles, the form of the Green 
functions simplifles at small scales such that the flow equations can be solved analyt- 
ically in this regime. The convergence of the Green functions towards these simple 
structures (as the factorized form ( p.l3|) ) in the infrared represents a new form of 



universality. 

We demonstrated the practical feasability of the method in the framework of 
the Wick-Cutkosky model, which represents a non-trivial bound-state problem with 
some known results derived via the Bethe-Salpeter equation. We checked the pos- 
sibility to continue the cm. energy s towards negative values (into the Minkowski 
regime), and the possibility to obtain both the mass of the lowest-lying bound state 
as well as the corresponding wave function with the help of the proposed combination 
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of numerical and analytic methods. 

The general applicability of these methods allowed us to investigate very different 
interactions such as QCD-motivated confining kernels. The relation between the 
Laplace transform of such kernels and confining nonrelativistic potentials including 
a constant term is actually an interesting subject for its own, we just briefly presented 
some general formulas on this relation. Again our results look quite promising, but 
this field deserves further studies. 

Clearly our approach can be improved and extended in a straightforward manner: 
The dependence of the four point function on the momentum variables wi,W2 can 
be restored, and fermions can be included [|l^. Furthermore the flow equation for 
the two point function or the exact propagator can be taken into account, and 
it is possible to go beyond the ladder approximation by including the six or even 
higher point functions. Also the consideration of gauge interactions will impose no 
intractable problems [0 : Of course the momentum cutoff has to be covariantized, 
which is straightforward since the UV cutoff A remains fixed throughout the whole 
procedure. It will turn out, however, that an additional one-loop counterterm is 
required |T^; this field is the subject of present investigations. 

In conclusion, besides a derivation and discussion of the flow equations we have 
shown that the bound-state problem is another field of practical applicability of this 
approach. It seems to be general and flexible enough to allow for studies of a wide 
range of phenomena in quantum field theory. 
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Figure Captions 

Fig. 1: Diagrammatic representation of the flow equations for the two and four 
point functions as in eq. The lines denote the full propagators F^^, 

encircled numbers the corresponding amputated point function, and the 
cross in a box an insertion of dkR^- 

Fig. 2: Plot of of the lowest-lying bound state in the Wick-Cutkosky model 
versus the coupling A. The short line denotes the known analytic result in the 
weak coupling limit (see eq. (|4.3|)), and the crosses known results for A = 1 
(eq. ( |4.5| )) and A = 27r (eq. ( |4.4| )). Our results are shown as error bars. 
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Fig. 3: Plot of M^, in the same model and for A = 1, versus the lattice size. The 
known result (eq. ( [4.5|) ) is indicated as a cross. 

Fig. 4: Plot of the logarithm of the Laplace transform of the wave function /(/) 
of the Wick-Cutkosky model in the weak coupling limit. The known result 
(eq. {\i.7f) ) is indicated as a straight line, and our results are normalized to the 
known one at / = 1. 
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Table 1 



m[GeV] 


M[GeV] 


1.5 


2.99-3 


2.5 


4.92-4.99 


3.5 


6.81-6.96 


4.5 


8.63-8.94 


5.5 


10.11-10.83 



Result for the bound state mass M for an interaction corresponding to a nonrel- 
ativistic potential ( |4.11|) and parameters as indicated below eq. ( [4.12| ). m denotes 
the different masses of the scalar constituents. 
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